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ABSTRACT 

We present a stacked weak-lensing analysis of a complete and volume-limited sample of 50 massive 
galaxy clusters at 0.15 < z < 0.3, based on deep, wide-field, high-resolution observations with Suprimc- 
Cam on the 8.2-m Subaru Telescope. Our new approach to selecting lensed background galaxies 
measures explicitly contamination by unlcnscd faint cluster and foreground galaxies, enabling us to 
define a secure sample of background galaxies lying red-ward of the cluster red sequence. Wc achieve 
a contamination fraction of just 1% in our background sample, from which we detect the tangential 
shear signal in the stacked data at a total signal-to-noise ratio of S/N = 32.7. The mass distribution 
reconstructed from the stacked data is nearly circularly-symmetric, and the tangential shear profile is 
remarkably smooth, justifying a one-dimensional radial analysis of the cluster mass distribution. The 
universal Navarro-Frenk- White model gives an excellent fit to the data, returning tight constraints 
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new results confirm and strengthen our previous claim that massive galaxy clusters have shallow 
density profiles that are broadly in line with predictions from numerical simulations. The sub- 10% 
statistical precision of our measurements is comparable with the differences between the predictions 
from numerical simulations in the literature. 

Subject headings: galaxies: clusters: general — gravitational lensing: weak — cosmology: observations 



1. INTRODUCTION 

In the cold dark matter structure formation paradigm 
galaxy clusters are the most massive, and most recently 
formed collapsed structures. A clear prediction derived 
from numerical simulations is that the density profile of 
galaxy clusters should therefore be less concentrated than 
that of less massi ve objects - galaxy groups, and in divid- 
ual galaxies (e.g IBullock et ajj|200lt IDolag et al.l I2004t 
INeto et al.l[2007tlDuffv et alJl2008t IZhao et alJ 12009^ 7" 

Gravitational lensing is a powerful probe of cos- 
mic matter density because the lensing signal is in- 
sensitive to the physical nature and state of mat- 
ter. The density profile of strong-lensing clusters 
is claimed to be more concentrated than predicted 
(iGavazzi et al.ll2003b iKneib et al.ll2003HBroadhurst et al.l 
l2008t lOguri et all 20091). ev en after correcting selection 
biases ([llennawi et al.l 12007). Interpretation of these re- 
sults is complicated by the small size of the observed 
samples, and the selection of the most spectacular strong- 
lensing clusters. On the theoretical side, the interpre- 
tation is complicated by the intrinsic scatter in profile 
shapes, and the rarity of massive cluster-scale dark mat- 
ter halos within the simulated volumes. 

Our weak-lensing study of a representative sample of 
30 X-ray selected clusters showed that their density pro- 
files agree better with predictions than those of spec- 
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tacular stro ng-lensing clusters (lOkabe et aLll2010l here- 
after OklO). lOguri et al.l (|2012l ) have also shown that less 
extreme strong-lensing clusters are in better agreement 
with predictions. Both studies employed stacking tech- 
niques to derive the mean properties of clusters, thus re- 
ducing the intrinsic shape noise of weak-lensing measure- 
ments, and eliminating biases due to the orientation and 
internal structure of clusters a nd the projection effe ct of 
large-scale structures (see also Uohnston et al"1l2007|) . 

In this paper, we stack a complete, volume-limited 
sample of 50 galaxy clusters that we have observed 
with the Subaru Telescope as part of the Local Clus- 
ter Substructure Survey (LoCuSS 5 ). The superb image 
quality, large field of vi ew, and large collect ing area 
of Subaru/Suprime-Cam (|Mivazaki et al.ir2002f) are ide- 
ally suited to accurate measurements of faint galaxy 
shapes th at are required for robust weak-lensing anal- 
ysis (e.g. lOklOf ). Our sample compris es all clusters 
from the ROSAT All Sky Survey c atalogs (jEbeling et al.l 
[19981 [20001 : iBohringer et al.|[200l that satisfy L x [0.1 - 
2.4keV]/£(z) 2 - 7 > 4.2 x 10 44 ergs _1 , 0.15 < z < 0.30, 
n H < 7 x 10 20 cm- 2 , and -25° < 8 < +65°. where 
E(z) = H(z)/Hq is the normalized Hubble expansion 

tion (jPopesso et al.ll2005h . 

In [J2]wc describe our data and analysis; in fj3]we ex- 
plain our results and compare with numerical simula- 
tions; in 2] we summarize our conclusions. We use the 
concordance ACDM model of Om,o = 0.27, Oa = 0.73 
and H = 100/t kms _1 Mpc _1 ijKomatsu et all 1201 If) . 
In this cosmology the virial over-density at the mean red- 
shift of our cluster sample, (z) = 0.23 is A v j r = 113.77. 
All error bars are 68% confidence intervals unless other- 
wise stated. 

5 http://www.sr.bham.ac.uk/locuss 



rate, and selecting on Lx I E{z) 2 ' 7 mimics a mass selec- 
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2. SUBARU DATA AND WEAK-LENSING ANALYSES 

Sensitive wide-field optical imaging of all 50 clus- 
ters is available from the Subaru/Suprime-Cam 67 . The 
observations, data reduction, and analysis will be de- 
scribed in detail in a future article. In brief, all 50 
clusters were observed through two filters, 46 through 
y/i'-band filters, 2 each through the V/Iq- and g/i'- 
band filters. Hereafter we refer to the bluer filter as V 
and the redder filter as i'. All data were reduced us- 
ing standard routines (|Okl0f ). The full-width half max- 
imum of point sources in the reduced V/i'-band frames 
is 0'.'6 < FWHM < 0"9 and 075 < FWHM < 077, respec- 
tively. Photometric calibration to < 10% precision in 
both filters was achieved via observations of Landolt 
standard stars, an d double checked agai nst SDSS/DR8 
stellar photometry (lEisenstein et al.ll201lT). We corr ected 
for galactic extinction (jSchlaflv fc Finkbc incr 201l|). 

We measure the sha pe of faint galaxies using a mod- 
ified version of lOklOf s pipelin e, which is based o n the 
IMCAT 8 implementation of the iKaiser et al.l (|1995l here- 
after, KSB) method. The main modification is that we 
calibrate the KSB isotropic correction factor for indi- 
vidual objects using a subs et of galaxies detect ed with 
high significance v > 30 ([Umetsu et al.l l2010f ) . This 
minimizes the shear calibration bias that is inherent to 
KSB+ methods in the pr esence of measurement errors 
(|Okura fc Futamasdl20ll . 

We confirmed the reliability of our shape measure- 
ments us ing simulate d data that were generated using 
GLAFIC ([Qguril 120101 ) with point spread functions de- 
scribed by the Moffat profile with a range of seeing 
(075 < FWHM < 171) and power indices (3 < /3 < 12), 
as described in lOguri et al.l ()2012f ) . We obtain a multi- 
plicative calib ration bias (m) and additive residual shear 
offset (c) (see iHevmans et al.l [2001 iMassev et~aD 120071 
for definitions) of \m\ 0.03 and |c| £ 2 x 10~ 4 , respec- 
tively, for FWHM ~ 077. 

We define a secure sample of background galaxies 
based on their color offset from the red sequence of early- 
type galaxies in the (V — color-magnitude plane 
for each cluster. In the absence of dust-reddening, in- 
trinsic scatter in galaxy colors, and measurement errors, 
selecting galaxies redder than the red sequence (here- 
after "red galaxies"; AC = (V — i 1 ) - (V - i') B so > 0) 
yields solely background galaxies. However, in reality a 
non-zero color cut is required to eliminate contamination 
by faint red cluster galaxies. In contrast, interpretation 
of "blue galaxy" (AC < 0) samples is complicated by 
star-formation. We discuss red and blue galaxies below, 

6 Based in part on observations obtained at the Subaru Obser- 
vatory under the Time Exchange program operated between the 
Gemini Observatory and the Subaru Observatory. The Gemini 
Observatory is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under a cooperative agreement with 
the NSF on behalf of the Gemini partnership: the National Sci- 
ence Foundation (United States), the Science and Technology Fa- 
cilities Council (United Kingdom), the National Research Coun- 
cil (Canada), CONIC YT (Chile), the Australian Research Council 
(Australia), Ministerio da Ciencia, Tecnologia e Inovacao (Brazil) 
and Ministerio de Ciencia, Tecnologia e Innovation Productiva (Ar- 
gentina) . 

7 Based in part on data collected at Subaru Telescope and ob- 
tained from the SMOKA, which is operated by the Astronomy 
Data Center, National Astronomical Observatory of Japan. 

8 http : / / www . if a. hawaii /kaiser / IM C AT 
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Fig. 1. — Top: Stacked reduced shear G+, as defined in §2, 
as a function of color offset from the stacked cluster red sequence 
(filled circles). Filled diamonds show the mean shear calculated 
after rotating all galaxies through 45° . Solid curves show the best- 
fit lensing kernel plus contamination model of the trend of G+ with 
AC The vertical dashed-dotted magenta line shows the color cut 
at which the fraction of contaminants is 1%. Middle: The color 
distribution of all galaxies at 21 < i' < 26 (black dashed curve). 
The width of the red sequence of bright (i 1 < 20) cluster members 
is shown as the green dotted curve, with the green vertical lines 
de-marking the la width of the red sequence. The magenta curves 
show the color distribution of contaminants in our model, and upon 
which the 1% contamination cut is based. Bottom: The mean 
distance with respect to the brightest cluster galaxy for 50 clusters, 
normalized by a uniform distribution. 

retaining only red galaxies for our main analysis and re- 
sults. 

The mean tangential distortion strength averaged over 
(i) all galaxies satisfying each color cut, (ii) cluster- 
centric radii of 0.1 /i _1 Mpc< r < 2.8 /i _1 Mpc, and (iii) 
all 50 clusters, G + = (((#+))), increases monotonically 
with AC for red galaxies (Fig. [TJ upper panel). We in- 
terpret the steep slope at AC ^ 0.3 as arising from con- 
tamination by faint red cluster members. Indeed, the 
mean cluster-centric radius of red galaxies is an increas- 
ing function of AC at small AC, supporting this inter- 
pretation (Fig.[l] lower panel). At AC ^ 0.4 we interpret 
the shallow slope of G+ as arising from a slowly increas- 
ing redshift of the faint red background population as 
AC grows. 

We therefore fit a two component model to the data. 
One component describes the color distribution of the 
contaminating population as a Gaussian centered at 
AC = 0, with standard deviation a. The other com- 
ponent describes the expected dependence of G+ on AC 
in the absence of contamination, and is based on the 
mean lensing kernel, D(AC) = (Di s /D s ), fo r galaxies in 
the C OSMOS photometric redshift catalog (jllbert et al.l 
12009? ) that matches each color cut. The model for 
the color-dependence of G+ is therefore: G+(AC) = 
A D(AC) (1 — Bf(AC)), where A is a normalization fac- 
tor (required to convert D into shear in a simple manner) , 
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B is the normalization of the Gaussian contaminant func- 
tion at AC = 0, and /(AC* > 0) = [1 - erf(AC/ v / 2cr)]/2. 
This model has three free parameters: A, B, and a. 

The best-fit model describes the data well (Fig. [TJ up- 
per panel). Interestingly, the best-fit standard deviation 
of the Gaussian is a = 0.21. This is larger than ex- 
pected based on the width of the cluster red sequence for 
bright galaxies and the photometric errors on faint galax- 
ies, suggesting intrinsic scatter in the color of faint red 
cluster galaxies. The model describes explicitly the de- 
pendence of contaminating fraction on AC. We conser- 
vatively adopt a limit of 1% on contaminating fraction, 
which translates into a color cut of AC > 0.475. We 
select galaxies redder than this cut for the results pre- 
sented in §3; the mean number density of these galax- 
ies is 5.3 ± l.Oarcmin -2 per cluster, where the uncer- 
tainty is the standard deviation among the 50 clusters. 
We therefore achieve a total stacked number density of 
266.3arcmin~ 2 . 

We applied the same methods to blue galaxies, de- 
scribing the contaminating fraction as /(AC < 0) = 
[1 + erf(AC/- v / 2cr)]/2. The model does not describe the 
blue galaxies well, likely due to the non-monotonic color- 
rcdshift relationship for these galaxies. We also note 
that the faint blue population appear to be preferen- 
tially found at large cluster centric radii (Fig. [1] lower 
panel), suggesting that blue galaxy contamination may 
be dominated by galaxies in the cluster outskirts. 

3. RESULTS 

We detect each of the 50 clusters at a typical peak 
signal-to-noise ratio of S/N ~ 4 in two-dimensional 
iKaiser k Squires! (|1993l) mass reconstructions. We also 
stack the shear catalogs in physical length units cen- 
tered on the respective BCGs and reconstruct the average 
mass distribution of the clusters. The mass density con- 
tours become progressively rounder and the significance 
of the weak-lensing detection increases as more clusters 
are stacked. The averaged mass map for the full sam- 
ple is remarkably symmetrical, with a pronounced mass 
peak of S/N -27.5 (Fig. [2]). 

Similarly, the stacked tangential shear profile achieves 
a high statisti cal precision (Fig. Our stacking pro- 
cedure follows lUmetsu et al.l (| 2 1 If) . In brief, we stack 
in physical length units, centering the catalogs on the 
respective BCGs, and stack across the radial range 
100 /i _1 kpc< r < 2.8 /i" 1 Mpc in 14 log-spaced bins. 
We detect the signal at S/N = 32.7 using the full covari- 
ance matrix to take into account projected uncorrelated 
large-scal e stru c ture and intrinsic elli p ticity noise (e.g. 
Hoekstral 120031: iHoekstra et al.l 120111: lOguri k Takadal 



201lUUmetsu et al.H2011t lOguri et aTll2012l ). computing 



the cosmic-shear co ntribution using th e non-linear mat- 
ter power spectrum (jSmith et al.ll2003t) for the WMAP7 
cosmology and the shape noise from the diagonal ma- 
trix. The 45°-rotated distortion component is consistent 
with a null signal, confirming that residual systematic 
errors are at least an order of magnitude smaller than 
the measured lensing signal. 

The observed density profile of massive galaxy clusters 
at z = 0.23 is clearly curved (FigEH) and is well-described 
by the parametric form proposed bv lNavarro et al.l ([19971 
so-called NFW profile): p cx a; _1 (l + x)~ 2 , where x = 



r/r s , and r s is the radius at which d\og p/dlogr = —2. 
We express our model fits in terms of the virial mass 
M V h= (47r/3) j o cr A v i r r v ; r 3 , and the concentration param- 
eter c v h= r V i r /r s , where A v ; r is the virial over-density 
and p cr is the critical density. We measure both param- 
eters to sub-10% precision (Tabled]), obtaining a best fit 
concentration parameter of c V i r = 5.4 ± 0.5. This result 
changes by just Ac V i r — 0.1 when we vary the number of 
bins between 8 and 18, change the inner radial cut from 
80 to 200 h~ 1 kpc or the outer radial cut between 2.5 and 
3.5 h~ 1 Mpc. The stability of our results under variations 
of the inner radial cut underlines the robustness of our 
new approach to selecting red galaxies, and the negligible 
level of (S x ) noted above. Finally, we note that our con- 



straints on concentration arc stable to Ac v 



0.2 with 



respect to increasing the color cut beyond AC > 0.475, 
and to fitting only to galaxies brighter than i' = 25. The 
constraints on M v ; r are stable to a few per cent under 
the same tests. 

We also investigate whether our stacking procedure 
introduces systematic errors. We use analytic NFW 
halos, that match the mass-concentration relation pre- 
dicted from numerical simulations, to create synthetic 
weak shear catalogs. These synthetic catalogs match 
the number density and the field of view of our Sub- 
aru data. We draw 300 samples of 50 clusters from the 
predicted cluster distribution, and stack the respective 
shear profiles in both physical length units and length 
units scaled to V2oo of each analytic NFW halo. We do 
not detect any bias in the measured mean concentration 
of the stacked clusters, obtaining (c/c tru th) = 1.02±0.07, 
and (c/ctruth) = 1-08 ± 0.07 for stacking in physical 
length units and re-scaled length units respectively. In 
both cases we obtain (Af/M tru th) = 0.96 ± 0.06, and the 
quoted uncertainties are the standard deviation on 300 
samples. The non-detection of a systematic error arising 
from stacking in physical units is consistent with OklO's 
result that their mass-concentration relations from indi- 
vidual and stacked clusters (using physical length units) 
are self-consistent. We also note that stacking in rcscaled 
length units weights the contribution of each cluster to 
each bin in a nonlinear and model-dependent manner: 

w cx 8A8 cx r| 00 cx M 2 qq . The impact of this weighting 
scheme deserves further detailed investigation in future 
studies. 

The test described above is idealized in that real clus- 
ters are aspherical, are embedded in the large-scale- 
structire, and contain baryons. Tests based on numer- 
ically simulated clusters should allow these issues to be 
investigated in detail in a future article. Here we con- 
duct a preliminary test using clusters extracted from 
the new "Cosmo-OWLS" si mulation that im p lemen ts the 
AGN model described in iMcCarthv et all (|2011[ ) in a 
400 /i _1 Mpc box, the weak-lensing catalogs were con- 
structed following iBahe et"al~l (|2012l ). The results arc 
consistent with the analytic NFW tests described above 
- the systematic error on the measurement of concentra- 
tion based on stacking in physical length units is sub- 
dominant to the statistical errors, and the central value 
on the estimated systematic error on stacking in rescaled 
length units is a factor of several larger than that on 
stacking in physical length units. 

Our statistical errors on concentration are com- 



4 




-500 



x [kpcft. ] 



x [kpc/i ] 



Fig. 2. — The projected mass distribution reconstructed from our weak-lensing catalogs: one typical cluster (ABELL0141; upper-left 
panel), 5 clusters (upper-right), 20 clusters (lower-left) and all 50 clusters (lower-right). Each panel spans 1.2 X 1.2 h~ 2 Mpc 2 , with contours 
starting at S/N = 3, and spaced at AS/TV = 2. A Gaussian smoothing scale of FWHM = 2arcmin is used in all panels (hatched region at 
lower right). 



TABLE 1 
Density Profile Models 



Model 


Shape parameter a 


A/vir 

(10 14 M Q /fc) 


c • b 


Af200 

(10 14 M G /fc) 


C200 b 




NFW 

gNFW 

Einasto 


7=1 

_ -, 27 +0.24 
7 ~ L " £l -0.37 

„ _ n 18S+ ' 062 
a - U.l»8_ . 058 


7 19 +u.m 

7 50+ ' 74 
'• OU -0.65 
7 40+0.86 
' -^ a -0.73 


5 41+ u - 4y 

A SS+0'86 
-0.80 


5 - 98 -o!38 

D - lo _0.44 
fi 15+0' 50 


4 22+ U ' 4U 
^•^ -0.36 
o 7Q +0.69 
°- ' 3 -0.69 
o on + 0.48 
°-° -0.66 


7.2/12 
6.6/11 
6.6/11 



a Parameter describing the shape of the mass density profile on small scales. 

b NFW- like concentration parameter denned by c^ FW = r&/r s , c gI ^ FW = (r&/r s )/(2 

and cf nasto = r A /r_ 2 . 
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parable with the differences between the predictions 
from different numerical simulations (Fig. |4}. In- 
deed the observed concentration parameter exceeds 
the predicted concentration of simulated clusters 
(iDuffy et alj|2008t iZhao et al.ll2009l: iBhattacharva et al.l 
120111: be Boni et al.ll2012D . iPrada et al.l (|2012T > recently 
noted that their mass-concentration relation from four 
simulations each comprising 10 7 particles agrees better 



with observations than previous work. Our best-fit c on- 
centration is slightly lower than that of IPrada et al.l . 

The NFW model fit (x^ in /d.o.f. = 7.2/12) does not 
motivate fitting more flexible models to our data (Ta- 
ble [I}. For completeness, we fit the generalized NFW 
(gNFW) and Einasto (1965) profiles. The former adds 
a free parameter 7 to the NFW profile: p oc x~ 7 (l + 
the latter describes the shape of the profile slope: 
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Fig. 3. — Stacked tangential shear profile of all 50 clusters in 
units of projected mass density, where different cluster and back- 
gro und galaxy redshifts galaxies are weighted by the lensing ker- 
nel fMa ndclbauin ct al. 2006; Okabe et aI.II20ld ; IQguri fc Takadal 
120111 ; lUmetsu et al .1120111 ). The projected radius is computed from 
the weighted mean cluster redshift (z c i uater ~ 0.23). The solid, 
dashed, dotted and dashed-dotted curves are the best-fit Navarro- 
Frcnk- White (NFW), singular isothermal (SIS), generalized NFW 
(gNFW) and Einasto profiles, respectively. The lower panel shows 
the result of the 45° test for systematic errors. 

dlogp/dlogr = — 2(?'/r_2) Q ; where r_2 is the radius at 
which the slope is isothermal. 

The best-fit gNFW profile is consistent with NFW, 
with 7 = 1.27^0 37. The best-fit Einasto profile has a = 
0.19 ± 0.06, consiste nt with numerica l simulations, e.g. 
a) = 0.175 ± 0.046 (IGao et all 120121) . and (a) = 0.183 
Navarro et all 120041 ). We also measure the inner slope 
of the best fit density profile models directly, obtaining 
/3(r = 0.01r 20u ) = -dlogp/dlogr = 1.1 for the gNFW 
and Einasto m odels, in good ag reement with ((3) ~ 1.1 
(jNavarro et al.ll2004t IGao et al.l l2012t ). This underlines 
that a constant power law slope is strongly disfavored; 
we obtain x m i n /d.o.f. = 44/13 for the singular isother- 
mal sphere (SIS) model [fi = 2). 

We examine the possible impact of adia batic contrac- 
tion o n the total measured density profile (|Gnedin et al.l 
120041 ) by introducing a central point mass contribution as 
an additional free parameter. We obtain an upper limit 
on the point mass of Mp j nt ^ 12 x 10 12 H^Mq, which is 
degenerate with the structural parameters of the smooth 
component in all models (NFW, gNFW, and Einasto). 
The best-fit mass and concentration parameters do not 
change significantly from those listed in Table 1 . The ex- 
cellent fit of the NFW model - that is based on numerical 
dark matter only simulations - to our weak-lensing data, 
and the results of adding baryons to the model (albeit in 
a simplified form) suggest the dark matter may not suffer 
adiabatic contraction by baryons in the cluster core. We 
will return to this topic in future more detailed study, 
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Fig. 4. — Stacked weak-lensing constraints on the mass and con- 
centration of a complete volume-limited sample of 50 galaxy clus- 
ters at (z) = 0.23. The white cross denotes the best-fit parameters 
and the contours show the 68.3%, 95.4%, and 99.7% confidence 
levels. Note that the predicted relations have all been converted to 
be consistent with our analysis. 

combining both strong- and weak-lensing constraints. 

We also checked whether the results are affected by 
adopting the BCG as the center of each cluster, by 
adding a n off-centering parame ter <jR a{{ to the models 
following IJohnston et all (|2007l ). The best- fit M vir and 
c v ; r arc unchanged, and we obtain an upper limit of 
cTR off < 29 /i _1 kpc. 

4. SUMMARY 

We have used sensitive high resolution observations 
with Subaru to measure the average density profile of 
a statistically complete volume-limited sample of 50 X- 
ray luminous galaxy clusters at 0.15 < z < 0.3. Care- 
ful treatment of systematic errors indicate that they 
are all smaller than the statistical errors. In particu- 
lar, we achieve just 1% contamination of the background 
galaxy sample by foreground and cluster galaxies, tests 
on simulated data indicate that our shape measurement 
multiplicative systematic error is m 0.03, and errors 
from choice of binning scheme are just a few per cent. 
When the signal from all 50 clusters is combined together 
we achieve a number density of background galaxies of 
266.3arcmin~ 2 ; ~ 200,000 galaxies in total. 

The shape of the stacked density profile is consis- 
tent with numerical simulations across the radial range 
100 /i _1 kpc —2.8 /i _1 Mpc. Specifically, we find no sta- 
tistical evidence for departures from the NFW profile. 
We constrain the mean mass and concentration of the 
clusters to sub-10% precision, obtaining c V i r = 5.41^°'^. 
This level of precision is comparable with the differences 
between the concentrations predicted by different numer- 
ical simulations, and therefore opens the possibility of 
discriminating between different simulations using obser- 
vational data in the near future. 

Our results emphasise the power of stacked weak- 
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lcnsing for constraining to high precision the average 
mass and shape of galaxy clusters. Imminent surveys 
such as Hyper Suprime-Cam on Subaru, the Dark En- 
ergy Survey, and KIDS, all hold much promise for simi- 
lar studies of less massive clusters, and clusters at higher 
rcdshifts. However significant advances on the precision 
that we have achieved here on massive low redshift clus- 
ters await future facilities such as LSST and Euclid. This 
is because only they will probe sufficiently deep and wide 
to improve on the number density of background galaxies 
on our sample of rare massive low redshift clusters, that 
is drawn from the ROSAT All Sky Survey. 
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